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The near zone energy-momentum tensor of a supersonic heavy quark jet moving through a 
strongly-coupled Af = 4 SYM plasma is analyzed in terms of first-order Navier-Stokes hydrody- 
namics. It is shown that the hydrodynamical description of the near quark region worsens with 
increasing quark velocities. For realistic quark velocities, v — 0.99, the non-hydrodynamical region 
is located at a narrow band surrounding the quark with a width of approximately 3/nT in the 
direction parallel to the quark's motion and with a length of roughly 10/-7rT in the perpendicular 
direction. Our results can be interpreted as an indication of the presence of coherent Yang-Mills 
fields where deviation from hydrodynamics is at its maximum. In the region where hydrodynamics 
does provide a good description of the system's dynamics, the flow velocity is so small that all 
the nonlinear terms can be dropped. Our results, which are compatible with the thermalization 
timescales extracted from elliptic flow measurements, suggest that if AdS/CFT provides a good 
description of the RHIC system, the bulk of the quenched jet energy has more than enough time 
to locally thermalize and become encoded in the collective flow. The resulting flow pattern close to 
the quark, however, is shown to be considerably different than the superposition of Mach cones and 
diffusion wakes observed at large distances. 



I. INTRODUCTION 

One of the most prominent experimental discoveries 
made at the Relativistic Heavy Ion Collider (RHIC) 
has been the suppression of highly energetic particles 
[TJ 13 [3J H], which suggests that the matter created at 
RHIC is a color-opaque, high density medium of colored 
particles where fast partons quickly lose energy by gluon 
emission [5J [5J [TJ [H|. Measurements of anisotropies in 
soft particle momentum distributions [TJ [2j [3j 0] have 
further indicated that soft degrees of freedom are ap- 
proximately thermalized. The degree of thermalization 
has been found to be considerably above the predictions 
obtained within perturbative quantum chromodynamics 
(QCD) [9| and, in fact, it seems to be compatible with 
the "perfect fluid" scenario where the strongly-coupled 
quark-gluon plasma (sQGP) has almost zero viscosity 
[10l [HJ [121 H31 Ell US]. This indicates that perturbative 
methods may fail at describing the microscopic proper- 
ties of the matter created in ultrarelativistic heavy ion 
collisions. 

An alternative approach for understanding jet suppres- 
sion in the sQGP relies on analyzing the similar prob- 
lem of a heavy quark moving through an AT — 4 su- 
per Yang-Mills (SYM) plasma. There the Anti-de Sit- 
ter/ Conformal Field Theory (AdS/CFT) correspondence 
[16j provides a way to compute observables at strong cou- 
pling on the CFT side, such as the stress tensor of the 
Af = 4 SYM plasma, starting from gravitational calcula- 
tions in five-dimensional AdS space. 

It is important to remark that one has to be careful 
when comparing QCD and Af = 4 SYM given that both 



the symmetries and the microscopic degrees of freedom 
are very different. The hope is that some of the differ- 
ences become unimportant in the strongly coupled limit 
where the system is to a good approximation a locally 
thermalized liquid. 

A notable result obtained using the AdS/CFT corre- 
spondence [TTJ was that the ratio between the viscosity 
r\ and entropy density s of an Af — 4 SYM plasma in 
the limit of large number of colors N c and large 't Hooft 
coupling A converges to a "universal limit" 



which has been conjectured [TB] to be the lowest value for 
this ratio obtainable in a physical system. Remarkably, 
hydrodynamic models find that the viscosity of the sys- 
tem created at RHIC is compatible, or even lower, than 
this limit [TTJEEIIIS]- 

The apparent local thermalization at RHIC suggests 
the possibility that the energy lost in jets will itself be- 
come locally thermalized, i.e., encoded in the local hy- 
drodynamic flow. In the immediate vicinity of the jet 
the variation of the thermodynamic quantities is prob- 
ably very large, which means that a hydrodynamic de- 
scription of the system is bound to fail in this region. 
The typical length scale associated with this region is, 
for instance, given by the sound attenuation length T = 
Arj/3uj < 1/ttT, where u> is the enthalpy of the medium 
(see [JS]). For distances scales larger than T, one ex- 
pects that the system can be described in terms of non- 
linear viscous hydrodynamics. Sufficiently far from the 
jet, however, the flow velocity and the gradients of the 
thermodynamic quantities are so small that linearized 
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hydrodynamic equations should provide an accurate de- 
scription of the system's dynamics. In the linearized 
limit, for certain energy deposition patterns, this gives 
rise to Mach cones, which can be readily connected to 
the fluid's equation of state [20 . In the absence of these 
conditions, the flow pattern is likely to be considerably 
more complex but will nevertheless result in a correlation 
between the flow and the jet. The possibility of detect- 
ing these effects in heavy ion collisions has recently been 
subject to intense theoretical and experimental interest 

[2]] |32l |23 121 12!3 123 



The degree of thermalization of the energy of the sup- 
pressed jets is a fundamental probe of the microscopic 
properties of the matter created at RHIC. Because of the 
apparent strong coupling of the system, non-perturbative 
methods such as AdS/CFT seem to be necessary. Re- 
cently, this method has been used to obtain the energy- 
momentum tensor of an infinitely heavy quark moving 
at constant velocity through a strongly coupled plasma 
[2S1 125] . The results have been analyzed in the far field 
limit [2911301 El E2J [331 IS] in terms of linearized Navier- 
Stokes hydrodynamics and, encouragingly, it has been 
shown that structures resembling Mach cones and diffu- 
sion wakes, with an angle comparable to that expected 
from the equation of state, do form in the far-field limit. 
In fact, in [M] it was shown that linearized Navier-Stokes 
hydrodynamics provides a very accurate description of 
the system at distances > 5/ ttT (in the limit of the large 
N c and A). 



As mentioned above, close to the quark linearized hy- 
drodynamics, and eventually the hydrodynamic descrip- 
tion, should fail to account for the energy- momentum 
tensor of the system. Distinguishing between non-linear, 
turbulent hydrodynamics, and the failure of hydrody- 
namics is impossible within the linearized regime. In this 
work we focus on the region very close to the quark and 
perform an analysis of the near quark energy- momentum 
tensor computed in [35] (see also [35] and [37]) using vis- 
cous, non-linear (first-order) Navier-Stokes hydrodynam- 
ics. We construct the flow vector by boosting the energy- 
momentum tensor in [35J to the local rest frame of the 
flow. This flow vector is used to construct the energy- 
momentum tensor corresponding to a viscous fluid having 
an equation of state and transport coefficients identical 
to those of strongly coupled Af — 4 SYM matter with 
the same flow vector. We then quantitatively examine 
the discrepancies between the two tensors and plot them 
as a function of the distance from the heavy quark. 



II. THE NEAR-QUARK ENERGY-MOMENTUM 
TENSOR 

The energy-momentum tensor of a heavy quark passing 
through an M = 4 SYM plasma at finite temperature T 
can be computed by considering metric fluctuations due 
to a string that is hanging down from the boundary of an 
AdS Schwarzschild (AdS-SS) background geometry [25] 
1 . In fact, in this limit the total action that describes the 
supergravity approximation to type IIB string theory in 
an AdS-SS background and the classical string is given 
by the sum of the following partial actions 

and 

Ang = -t^—, [ \f-G { °}d a Xnd X»d 2 <T, (3) 
Zna J " 

where L is the radius of AdS 5 , G 5 = ttL 2 /2N 2 , a' — 

L 2 /V% Giiv is the total metric, and G®} is the metric of 
the unperturbed AdS-SS black hole (without effects from 
the string), which can be obtained from 

d*=*(-g(,)dl* + d3* + ^) (4) 

where z goes from at the AdS boundary to zq — 1 /irT 
at the black hole horizon (T is the Hawking temperature 
associated with the black hole) and g(z) = 1 — (z/zo) 4 . 
The string coordinates X^(a, r) in the Nambu-Goto ac- 
tion in Eq. ([3| are chosen in such a way that the string 
endpoint (which corresponds to the heavy quark in the 
4-dimensional boundary) moves at constant speed v and 
no energy flows from the horizon into the string |38l 139] . 

Minimizing the action S with respect to G leads to 
the full set of Einstein's equations. It is sufficient for 
our purposes here to consider instead the linearized Ein- 
stein's equations for the metric fluctuations h^ u , which 

are defined via G M „ = G^v+h^. It can be shown [251 [35] 
that the contribution from the moving quark to the total 

energy- momentum tensor is T quar ^ = ^yj i^ v s Qi where 
the tensor Q is obtained by expanding h in powers of 
z near the boundary, i.e., h ~ Qz 4 . The total energy- 
momentum tensor in the lab frame that describes the 
near-quark dynamics of the plasma was computed by 
Yarom [35] and it reads 

T* = P diag{3, 1, 1, r 2 } + AT^(x u r). (5) 



Throughout this paper the units are H = c = ks — 
1 and 4-vectors are denoted by capital letters, e.g., 
jjiJ. — (jjQ.ijy Also, the Minkowski metric g^ v = 
diag(— , +, +, +) is used. 



Fluctuations of the string can be neglected when M/T 2> vA, 
where M is the quark's mass |23l 1381 140| . In this paper we con- 
sider only this classical limit. Assuming that T = 0.5 GeV and 
A = 3ir, one sees that the classical string dynamics is a much 
better approximation for bottom quarks than for charm quarks. 
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where the explicit form of AT MZ , is 

v (r 2 (-5 + I3v 2 - 8v 4 ) + (-5 + nv 2 )xf) x x 



AT tt = a- 



72 (r 2 (l-v 2 ) + z 2 ) 5/2 



AT tXl = -a 



v 2 (2xj + (1 - v 2 )r 2 ) xi 



24(r 2 (l-v 2 ) + x 2 ) 



2^/2 



AT,, 



AT XlXl = a- 



(1 - v 2 )v 2 (llxj + 8r 2 (l - v 2 )) r 
72(r 2 (l -v 2 ) + x 2 f' 2 

r 2 (8 - 13v 2 + 5t> 4 ) + (11 - bv 2 )x\) x x 
72 (r 2 (l - v 2 ) + x 2 f' 2 



AT Xir = a 



w(l-t; 2 )(8r 2 (l- U 2 ) + lla; 2 )r 



72 (r 2 (l-w 2 ) + x 2 ) 



2^5/2 



AT, 



v(l - v 2 ) (5r 2 (l - w 2 ) + 8xf) xi 

t 

72 (r 2 (l - v 2 ) + x\f /2 



ATaa = 



r 2 v(l — v 2 )x\ 



9(r 2 (l-v 2 ) + x() 



(0) 



and a = ^qvX^T 4 . In the equations above we used the 
dimensionless cylindrical coordinates xi — X\ ttT (where 
Xi = X — vt is the comoving coordinate of the quark), 
r = X p TrT, and also defined 7 g = — v 2 . Moreover, 

P = (N 2 - 1) vr 2 T 4 /8 + 0(N°) is the pressure of the 
ideal SYM plasma jH]. Note that the disturbances in 
the medium caused by the string are parametrically of 
0(T 2 /X 2 ) in the near-quark region [35], while the vac- 
uum Coulomb stress (not included here) is 0(1/ X 4 ) [21)] . 

It is assumed throughout the derivation of Eq. ([5| 
that the metric disturbances caused by the moving string 
are small in comparison to the AdS§ background metric. 
Therefore, this result is correct as long as this condition 
is fulfilled. In fact, since AT M „ scales inversely with the 
total distance from the quark, the region where the con- 
dition AT^/Pq < 1 (or, equivalently, h small in compar- 
ison to G(°>) holds can be taken to be arbitrarily small 
as long as the limit where N c — > oo and A is large is 
employed. However, in order to evaluate the relevance 
of this approach to heavy ion collisions, one may choose 
N c = 3, A = 3n (a s = 0.25), and -f q = 10. Using these 
parameters, one obtains that the magnitude of the en- 
ergy density disturbances caused by the string along the 
jet axis, AT t4 / (3P ) |r=o, is < 1 when > 1/(ttT) 

(which is comparable to uncertainty principle bound on 
the mean free path discussed in [3]). 



III. COMPARISON TO HYDRODYNAMICS 

We shall now compare the tensor in Eq. [5] to the one 
that represents a solution to the relativistic Navier-Stokes 
equations with an arbitrary flow. We underline that the 
quantitative values for the flow and thermodynamic quan- 
tities resulting from our analysis should be taken with 
caution. In fact, one would expect that the validity of 
a low-energy hydrodynamic description of the near-field 
region worsens as one becomes closer to the heavy quark. 
However, this analysis is necessary to ascertain the do- 
main of validity of (general) hydrodynamics, as opposed 
to non-equilibrium field configurations (up to the "op- 
posite extreme" of fully anisotropic classical solutions to 
the field equations). As we show in this paper and sub- 
sequent work |42j . ascertaining this domain of validity 
of hydrodynamics is crucial for the development of the 
phenomcnological properties of AdS/CFT based models. 

In hydrodynamics, the Knudsen number Kn is defined 
as the ratio between the mean free path Imfp and a 
characteristic spatial dimension of the system q. Hydro- 
dynamics is applicable when Kn = Imfp/q < 1. In 
conformal field theories at finite temperature, the only 
dimensionful parameter is given by the temperature T 
and, thus, both Imfp and q should be proportional to 
1 jT . However, the mean free path is not a well defined 
quantity in M = 4 SYM theories at very strong coupling. 
Nevertheless, one can still define an effective Knudsen 
number as 



K N =T 



V-5 



S 



(7) 



where T = l/3irT is the sound attenuation length (to 
leading order in N c ), Si = —T^ is the momentum den- 
sity, and S = Vd 2 ". Note that this quantity does not 
depend on N c and A but it strongly depends on v. It is il- 
lustrated in Fig. [T]how Km changes with v. For v — 0.75, 
Kn <C 1 practically everywhere but in a small region of 
radius ~ 1/nT that surrounds the quark. However, when 
v takes more realistic values, say, v = 0.99, > 1 in 
a narrow region located at —3/irT < X\ < 3/irT and 
X p > 0. Therefore, one can expect that in this region 
the system cannot be described by hydrodynamics. 

Another way to verify whether the system described 
by the energy-momentum tensor in Eq. ^ admits a hy- 
drodynamic description is to compute its flow vector , 
which can be obtained by boosting to the Landau frame 
where (T^)x = (denoted henceforward by brackets and 
subscript L, (...) L ). Note that, unless the system is a co- 
herent field where the phase velocity is equal to the speed 
of light (such would be the case of an electromagnetic 
wave), this transformation is always possible. It can be 
accomplished by solving a system of two equations for the 
two space-like components of ?7 M , U\ and U r (Ug = 0), 



IT, 



) L =A?T*AZ = 0, 



(8) 



where is a general coordinate dependent Lorentz 
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that Eq. pi becomes 
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FIG. 1: (Color online) Effective knudsen number for v = 0.75 
(left panel) and v = 0.99 (right panel). Note that different 
scales were used in the left and right panels. Ratio values 
outside of the limits were color-coded as the limits. 



transformation 



A = 



U 1 



u 1 

u<s>u 7 
u 2 



( 7 -l) 



(9) 



where 7 = U° — yl + U 2 . Using the representation 



T Y = 



e -S T 
-S f 



where e — Tq and ry = T? — 1,2,3), we obtain 



U = 



1 



7£ 



S T -U 



u®u 7 
u 2 



(7-1) 



(75 



-fU 



(in 

For finite N c and A this equation can only be solved nu- 
merically. However, for very large N c and large A one 
can use that, to leading order in 1/N C , 7 « 1, e as 3Po, 
and t rs Pq. Using these approximations, one can then 
obtain that 



U 



s 

4iV 



(12) 



In deriving the equation above we used that ^/X/N 2 <C 1. 
The fact that \U\ - 0(\f\/N 2 ) in the large N c , A limit 
implies that nonlinear terms in the hydrodynamic de- 
scription arc sublcading contributions that can be ne- 
glected. Thus, if the system's dynamics can be described 
by hydrodynamics, these equations have to be linear for 
the present string theory setup obtained in the large N c , 
A limit [33]. However, at finite N c and A nonlinear effects 
are expected to be relevant. These nonlinear effects can 
only be properly taken into account by incorporating sub- 
leading 1/N C corrections. Because these corrections have 
yet to be calculated, we shall extrapolate the validity of 
our leading order results and make a comparison between 
the near-quark energy-momentum tensor in Eq. (|5 1 and a 
hydrodynamic ansatz using the phcnomenologically rele- 
vant set of parameters iV c = 3 and A = 3tt. In this work, 
to estimate both the non-linearities and non-equilibrium 
effects in the near-zone, we shall use the non- linear ansatz 
in Eq. (tTTl) rather than Eq. (pi. 



The energy-momentum tensor of a viscous fluid con- 
taining no conserved charges can be parametrized using 
the energy density p, pressure p, and energy flow vector 



T, 



NS 



(p + p) U^U V +pg^ 



n„ 



(13) 



(in our metric convention U^U^ = — 1 ). n M!/ is the shear 
tensor, which contains the corrections to local thermal- 
ization that can be expanded in terms of the gradients 
of U^. To first-order in these gradients (Navier-Stokes 
limit), this tensor is given by [20] 

n"" = -^(d^tr + d u u» + u^u a d a u v 

+ U v U a d a U^) + \j 1 ^ {d a U a ), (14) 

where = g^ v + U^U V is the local spatial projector 
and 77 = (N 2 — l) 7rT 3 /8. Note that we kept the nonlin- 
ear terms in C/ M in the definition of H^. 

Once the solution to U has been found (by solving 
the transcendental equation (111 numerically), the shear 



(10) tensor n M „ and the energy-momentum tensor T^ s can 



be computed from Eqs. (14 1 and (13). The continuity of 
Ti^v confirms the stability of our numerical solution and 
the existence of only one physically relevant set of U. 
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For an inviscid fluid, the form of the transformed tensor 
will be 



K) £ =diag{p,p, p,p}. 



(15) 



If the first-order Navier-Stokes approximation provides a 
full description of the system, 



(T Y ) — (t ns ) — (T id ) 
rp/I V / l v ft") l 



inn! 



(16) 



where 



• (II'"') L is the shear tensor of Eq. (14 1 transformed 
to the Landau frame. 



p in Eq. ( 15 I is given by 



p{x 1 ,r) = (T y ) L -(n 00 ) L (17) 
• The equation of state is the CFT one, p = p/3. 



Note that Eq. ( 15 ) can in general be very different than 
the "background" tensor (T p i asma ) = P diag{3, 1, 1, r 2 }, 
since it is defined in the frame comoving with the flow 
generated by the jet. 

The deviation from Navier-Stokes hydrodynamics can 
then be quantitatively investigated by defining the tensor 



7 _ tt 

a tfC J ™ ix U!/ 



(18) 



and studying, in particular, discrepancies between 
(Z^) L and {Tjf v ) L . If {Tj v ) L corresponds to a solution 
of the Navier-Stokes equations, we should get 



and 



(Z n ) L = {Z 22 ) L = (Z 33 ) J 



IV. RESULTS 



(Z 00 ) L , (19) 



(20) 




FIG. 2: (Color online) Contour plot of the magnitude of the 
non-relativistic flow velocity V for v = 0.99. 



u 




In this preliminary analysis, we will concentrate on the 
case where A = 3n, N c =3, and the quark velocity in 
the lab frame is v — 0.99. We would like to remark at 
this point that 1 / ttT is the natural uncertainty principle 
bound on spatial resolution at temperature T in an ul- 
trarelativistic plasma and provides a lower bound on the 
mean free path and hence viscosity [5]. Within a radius 
1 /ttT of the quark jet our semiclassical dissipative fluid 
analysis is expected to fail and non-equilibrium and field 
coherence effects may dominate the dynamics. 

The plot seen in Fig. [2] immediately illustrates that 
the flow nearby the quark is very different than a super- 
position of a Mach cone and a diffusion wake (V is such 
that U = jV). Note that V 2 < 1 everywhere. This 
structure is expected to be somewhat stable with respect 
to v since viscosity results in a length scale below which 
the fluid "sticks" to and co-moves with the quark. In 



Xi n T 

FIG. 3: (Color online) Contour plot of the collective flow U 
for v = 0.99. 



a steady state solution, even a small viscosity can pro- 
duce an observable disturbance preceding the quark. In 
the hydrodynamic simulations of [3D 133 HS] the fluid is 
ideal. Hence, the "shock" -like nature of the quark distur- 
bance is recovered and the energy density is qualitatively 
different. 

In Fig. [3] we superimposed the magnitude and the vec- 
tor plot of the collective flow U. As the quark moves 
forward, matter is being transported by collective flow in 
the direction opposite to quark motion. This effect is, in 
itself, not completely surprising since perturbative calcu- 
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FIG. 4: (Color online) Comparison between the full numer- 
ical solution for U and the leading order expression in Eq. 
(H, i.e., \U-S/4P \ for v = 0.99. 



lations |47j admit both momentum flow out of the quark 
(if repulsive channels dominate) and into the quark (if 
attractive channels dominate) . Furthermore, in a system 
with no conserved charges, collective flow and heat diffu- 
sion can compensate each other since the effect of both is 
moving energy-momentum throughout the system. This 
picture, is however very far from the point-like energy 
source considered in the "text-book" Mach cone example 

of mi- 

In Fig. [4] we compare the full numerical solution of 
Eq. Q and the leading order result in Eq. ( |T2| ). This 
plot indicates that the system's dynamics becomes linear 
outside the narrow region between — 1 /ttT < Xi < 1/ irT 
and < X p < 6/ttT. This is in agreement with the result 

in Fig. [^because the magnitude of V (or U) becomes very 
small outside this region. 

Finally, Fig. [5] shows the discrepancy between the 
Navier-Stokes tensor and T^ v in the Landau frame, which 



are parametrized by the tensor Z^ (Eq. 18). As can be 
seen, the magnitude and structure of the discrepancy re- 
mains somewhat the same across all components of Z^. 
Outside the region where —3/irT < Xi < 3/ttT and 
< X p < 10/ttT (the quark is at X x = X p = 0), the 
system can be described by the first-order Navier-Stokes 
energy-momentum tensor to an accuracy of roughly 90%. 
Inside that region, however, the discrepancy quickly rises 
and diverges. Note that in the region where Navier- 
Stokes provides a good description of the system the nu- 
merically calculated velocity flow U matches the leading 
order in Eq. (12). This implies that, in this region, the 



system is better described by the linearized version of the 
first-order Navier-Stokes equations. 

The detailed structure of the divergence is numerically 
difficult to examine but its sudden onset suggests the 



presence of coherent Yang-Mills fields, where deviation 
from hydrodynamics is at its maximum. Coherent fields 
do not of course have a comoving frame, but they can 
be transformed into a frame (which will in general be 
different than the frame set by as calculated here) 
where the energy-momentum tensor given by 



frpcoherent\ 



diag{p',p',0,0} 



(21) 



(where p' and p' are taken to mean energy and momen- 
tum flow density, not "true" thermal densities and pres- 
sures). Here, the X% direction is taken to be parallel to 
the Poynting vector. The Z u calculated for this field can 
go to zero, which results in divergences in ratios of Z's 
as seen in Fig. [5] This indicates that coherent Yang- 
Mills fields are present in the vicinity of the quark. In 
fact, note that Yarom's AT M „ (j6| does seem to have a 
Lorentz contracted Weizsacker- William-like form but its 
<D(T 2 /X 2 ) parametric dependence is quite distinct from 
the vacuum Coulomb behavior 1421. 



Our results for the components of Z do not signifi- 
cantly change if we increase the intensity of the coupling 
to A = 67T. However, at lower quark velocities the Navier- 
Stokes equations provide an accurate description of the 
system's dynamics at distances > 1.5/nT (see Fig. |6j. 
These results coincide remarkably well with the limit es- 
timated through the uncertainty principle [S]. 

Fig. [7] confirms that stays smooth and finite 

throughout, and the divergence is due to a maximal 
breakdown of isotropy (one of the Zu going to zero). 
Furthermore, Fig. [7] shows that in the regime in which 
hydrodynamics works, is a perturbatively small cor- 
rection to the energy- momentum tensor. The regime 
where Tl^ becomes significant approximately coincides 
with the regime in which the Navier-Stokes description 
stops to be a good description of the physical system. 
Thus, it is unlikely that higher gradient corrections to 
Navier-Stokes equations (such as [1H1 HHJ [50] ) will make 
a significant difference in our results. 

The contours of our figures look non-trivial and no 
easy interpretation of them can be given here, especially 
as they are highly susceptible to 1/N C corrections (as re- 
marked in the previous section). The fine detailed struc- 
ture in the figures is sensitive to the interpolation nu- 
merics and only the global features are numerically sta- 
ble. We remark, however, that the size of the "inner re- 
gion" is very similar in all the figures and also agrees well 
with the estimate for the "Knudsen boundary" in Fig. [T] 
It is therefore worth underlining how the results agree 
with the expectations described in the previous sections. 
The system is neatly defined by a Knudsen region where 
—3/ttT < X\ < 3/ttT. Inside this layer, non-equilibrium 
deviations from the hydrodynamic description dominate. 
On the other hand, outside this region nearly ideal and 
linearized hydrodynamics provides an acceptable descrip- 
tion of the system. 
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FIG. 5: (Color online) Tests of local equilibrium for the elements of the energy-momentum tensor in the Landau frame. 
is defined in Eq. [l8]and v = 0.99. Ratio values outside of the limits were color-coded as the limits. 



CONCLUSIONS 



Encouragingly, the scale at which the energy deposited 
by the jet thermalizes is comparable to the thermaliza- 
tion time constraints obtainable from anisotropic flow 
measurements which suggests that in the strongly 
coupled regime the nearly inviscid microscopic dynam- 
ics responsible for creating large V2 such as those seen 
at RHIC also has the potential to transform the energy 
of quenched hard jets into flow of soft particles. The 
short timescale seen in the present work indicates that, 
provided strongly coupled AdS/CFT is indeed qualita- 
tively similar to sQGP, the energy deposited by jets in 
the system becomes to a very good approximation locally 
thermalized. 



Second-order hydrodynamic corrections [45] may 
somewhat change the results presented here. However, 
as argued previously, the contribution of the relaxation 
time to collective dynamics might be not so large. Ana- 
lyzes of other systems within the AdS/CFT framework, 
such as |51j . confirm this. An analysis of the presently 
investigated solution within the second order framework, 
analogously to that presented in [51] 122] for a Boost- 
invariant fluid, is currently in progress. 

It should be noted that there is to date no universally 
agreed framework for relativistic hydrodynamics beyond 
leading order, with several higher order frameworks pos- 
sible (see, for example, [48, 49, 50 and references within 
these works) . This still unresolved ambiguity might hold 
the key to explaining some puzzles encountered in trying 
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FIG. 6: (Color online) Test of local equilibrium for the ele- 
ments of the energy-momentum tensor in the Landau frame 
for v — 0.75 and A = 3ir. Zoo is defined in Eq. |18| Ratio 
values outside of the limits were color-coded as the limits. 




X\ 7T T 

FIG. 7: (Color online) The relative importance of the viscous 
tensor with respect to the energy-momentum tensor in the 
Landau frame, yl^/i/ njL/ corrected for double- 

counting. Here v = 0.99. 

to model heavy ion collisions using non-ideal hydrody- 
namics (the numerical instabilities encountered in |12j . 
and the very low viscosity encountered in |15j). For a 
discussion of second-order viscous hydrodynamics in con- 
formal field theories at finite temperature see [53J [52] . 

We hope that AdS/CFT techniques might help in clar- 
ifying the situation by providing a bonafide "analytical" 



solution for a strongly coupled relativistic quantum sys- 
tem with none of the approximations, whose validity is 
questionable at strong coupling, that usually go into the 
derivation of transport equations. For such a clarifi- 
cation, an analysis of AdS/CFT inspired solutions us- 
ing methods similar to those described here is necessary. 
This effort could, independently from AdS/CFT's appli- 
cability to heavy ion collisions, lead to a better compre- 
hension of the general and fundamental problem of the 
dynamics of strongly-coupled relativistic systems. 

In conclusion, in this work we have compared the an- 
alytically obtained energy-momentum tensor of a very 
massive quark moving through a strongly coupled Af = 4 
SYM plasma with a non-ideal tensor given by the first- 
order Navier-Stokes ansatz. We have found that the en- 
ergy deposited by the quark is already thermalized in the 
region where \Xi\ > 3/irT and X. p > 1/ttT, which de- 
fines a length scale that is considerably shorter than the 
time-scale of the jet traveling throughout the medium, or 
the medium's hydrodynamic evolution. Thus, if strongly 
coupled Af = 4 SYM is an appropriate description of 
RHIC physics, we have every reason to expect that in the 
strongly coupled region the missing jet energy has been 
thermalized. The thermalization timescale is (indeed, the 
lifetime of the heavy ion system), slightly smaller than 
the timescale found in Thus, in view of the re- 

sults present here, one can say that linearized first-order 
Navier-Stokes hydrodynamics provides an accurate de- 
scription of the Af = 4 SYM plasma surrounding the 
heavy quark down to distances of about 3/ttT. A quan- 
titative hydrodynamic analysis of the energy deposited in 
jets at smaller distances needs to be treated nonlinearly, 
as in gSJUS]. 

We expect that further progress in examining strongly 
coupled Af — 4 SYM systems will better clarify the ther- 
malization timescale and dynamics of strongly coupled 
field theories. 
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